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(57) Abstract 

A system and method non-invasive biomedical optical imaging 
and spectroscopy with low-level light is described. The technique 
consists of a modulated light source (120) coupled to tissue 
(100) of a patient to introduce excitation light. Fluorescent 
light emitted in response to the excitation light is detected with 
sensor (148). The AC intensity and phase of the excitation and 
detected fluorescent light is provided to a processor (160) operative! y 
coupled to sensor (148). Processor (160) employs the measured 
re-emission kinetics of excitation and fluorescent light to "map" 
the spatial variation of one or more fluorescence characteristics of 
the tissue (100). The fluorescence characteristic may be provided 
by exogenous contract agents, endogenous fluorophores. or both. 
The variations is determined by solving frequency domain diffusion 
equations at a number of designated points in the tissue as part 
of a recursive estimation algorithm. Processor (160) generates 
an imaging signal in accordance with the spatial variation of the 
fluorescence characteristic for provision to an output device (164). 
The output device (164) displays an image corresponding spatial 
variation of the fluorescence characteristic which corresponds to 
tissue (100) to aid in the detection and diagnosis of disease. 
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FLUORESCENCE LIFETIME- BASED IMAGING AND SPECTROSCOPY 
IN TISSUES AND OTHER RANDOM MEDIA 

BACKGROUND OF THE INVENTION 

The present invention relates to spectroscopic imaging of heterogeneous light 
scattering media, and more particularly, but not exclusively, relates to in vivo imaging 
of biologic tissue by mapping a fluorescence characteristic of the tissue through the 
detection of light emitted in response to excitation light from a time-varying light 
source. 

The early detection of disease promises a greater efficacy for therapeutic 
intervention. In recent years, non-invasive techniques have been developed which have 
improved the ability to provide a reliable and early diagnosis of various afflictions by 
detecting biochemical changes in the tissue of a patient. For example, Magnetic 
Resonance Imaging (MRI) has successfully monitored the relaxation of spin states of 
paramagnetic nuclei in order to provide biomedical imaging and biochemical 
spectroscopy of tissues. Unfortunately, the complexity and expense of MRI diagnostics 
limit its availability - especially as a means of routine monitoring for disease. 

Another powerful analytical technique with an increasing number of 
applications in the biological sciences is fluorescence spectroscopy. These applications 
include biomedical diagnostics, genetic sequencing, and flow cytometry. To date, there 
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are several industrial and academic institutions developing fluorescent and 
phosphorescent compounds for observing pertinent metabolites and environmental 
conditions, such as Ca**, pH, glucose, p0 2 , and pC0 2 . With the development of dyes 
and photodynamic fluorescent agents which excite and re-emit in the near-infrared red 
(NIR) wavelength regime, non-invasive detection of diseased tissues located deep 
within tissues may also be possible since red excitation and re-emission light can travel 
significant distances to and from the tissue-air interface Wilson et aL Time- 
Pepgndent Optical Spectroscopy and Imaging for Bi o medical App lication*: 80 
Proceedings IEEE pp. 918-30 (1992)). 

As exemplified by U.S. Patent Nos. 5,421,337 to Richards-Kortum et al. and 
5,452,723 to Wu et al., several investigators have suggested various procedures to 
differentiate diseased and normal tissues based on fluorescence emissions through non- 
invasive external measurements or minimally invasive endoscopic measuring 
techniques. Unfortunately, these procedures generally fail to provide a viable spatial 
imaging procedure. One reason imaging based on fluorescence has remained elusive is 
that meaningful relational measurements of fluorescence characteristics from a random, 
multiply scattering media, such as tissue, are difficult to obtain. For example, 
fluorescent intensity, which is a function of the fluorescent compound (or fluorophore) 
concentration or "uptake," is one possible candidate for imaging; however, when this 
property is used in an optically dense medium, such as a particulate (cell) suspension, 
powder, or tissue, the local scattering and absorption properties confound measured 
fluorescent intensities. 

Besides intensity, other properties of selected fluorophores such as fluorescent 
quantum efficiency and lifetime are also sensitive to the local biochemical environment. 
As used herein, "fluorescent quantum efficiency" means the fractional number of 
fluorescent photons re-emitted for each excitation photon absorbed or the fraction of 
decay events which result in emission of a fluorescent photon. "Fluorescent lifetime," 
as used herein, is defined as the mean survival time of the activated fluorophore or the 
mean time between the absorption of an excitation photon and re-emission of a 
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fluorescent photon. Like intensity, measurement of these fluorescence characteristics 
is often limited to well-defined in vitro applications in the research laboratory or in 
flow cytometry where issues such as scattering, absorption, and changing fluorophore 
concentrations can be controlled or measured. Moreover, these limitations generally 
5 preclude meaningful fluorescence-based imaging of hidden tissue heterogeneities, such 
as tumors or other diseased tissue regions which cannot be detected by visual 
inspection. 

Thus, a need remains for a technique to non-invasively image multiply 
scattering tissue based on one or more fluorescence characteristics which does not 
10 require extensive information about intrinsic optical properties of the tissue, and takes 
advantage of the contrast capability offered by fluorescence yield and lifetime 
characteristics to aide in the identification of tissue heterogeneities. The present 
invention satisfies this need. 
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SUMMARY OF THE INVENTION 

The invention relates to spectroscopic imaging of heterogeneous, light 
5 scattering materials. Several aspects of the invention are novel, non-obvious, and 
provide various advantages. While the actual nature of the invention covered herein 
can only be determined with reference to the claims appended hereto, certain features 
which are characteristic of the present invention are described briefly as follows. 

One feature of the present invention is a technique for imaging a heterogeneous 

10 light scattering material. This process includes exposing the surface of a material to 
light from a light source and detecting an emission in response. A spatial variation of a 
fluorescence characteristic of the material is determined as a function of the emission 
with a processor. The spatial variation may be characterized by a set of values 
representative of the fluorescence characteristic as a function of position. An image is 

15 generated in accordance with the spatial variation which corresponds to the 

heterogeneous composition of the material. This technique may be applied in vivo to 
biologic tissue using external or endoscopic instrumentation to detect heterogeneities 
indicative of disease. The technique may include the introduction of a fluorescent 
contrast agent into the material. The fluorescence characteristic detected may be 

20 fluorescence lifetime, fluorescence quantum efficiency, a fluorophore absorption 
coefficient, fluorescent yield (a function of fluorescent quantum efficiency and 
fluorophore absorption), or another fluorescence characteristic known to those skilled 
in the art. 

In another feature of the present invention, the spatial variation of a light 
25 scattering material with a heterogeneous composition is determined by establishing an 
estimate of the optical property or fluorescence characteristic variation, determining a 
calculated emission from the material as a function of the estimate, and comparing the 
calculated emission to a detected emission to determine a corresponding error. The 
estimation of the variation is modified, the calculated emission re-determined with this 
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modified estimate, and the comparison repeated until the error reaches a desired 
minimum. An image of the material is generated from the modified estimate which 
corresponds to the heterogeneous composition. 

Accordingly, it is one object of the present invention to map a fluorescent 
property of a light scattering material which varies with the heterogeneous composition 
of the material to generate a corresponding image. 

It is another object of the present invention to provide a spectroscopic technique 
for non-invasively monitoring fluorescent properties of hidden tissue volumes in a 
living organism and to monitor selected metabolites of an organism in vivo. 

Yet another object is to provide a fluorescence imaging system and process to 
identify diseased tissue using endogenous or exogenous fluorophores as a contrast 
agent. This contrast may be provided by a concentration, lifetime, or quantum 
efficiency difference between normal and diseased tissues. 

It is still another object of the present invention to provide an imaging technique 
and algorithm which is based on contrast of an optical property that is independent of 
local fluorophore concentration. 

Further objects, features, aspects, and advantages of the present invention will 
become apparent from the drawings and description contained herein. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a schematic illustration of a system of one embodiment of the present 
invention. 

5 FIG. 2 is a flow chart of a process utilizing the system of FIG. 1 . 

FIG. 3 is a schematic representation of a tissue phantom arrangement used to 
demonstrate various aspects of the present invention. 

FIGS. 4-7 graphically depict selected properties of equations used in the present 
invention. 

10 FIGS. 8 and 9 graphically depict convergence of simulated determinations of 

the spatial variation of fluorescent yield and lifetime, respectively, utilizing one 
embodiment of the present invention. 

FIGS. 10-14 are images obtained from experimental examples 1-3 of the present 
invention. 

15 FIG. 15 is a schematic illustration of a system of an alternative embodiment of 

the present invention. 
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DESCRIPTION OF PREFERRED EMBODIMENTS 

For the purposes of promoting an understanding of the principles of the invention, 
reference will now be made to the embodiment illustrated in the drawings and specific 
5 language will be used to describe the same. It will nevertheless be understood that no 
limitation of the scope of the invention is thereby intended. Any alterations and further 
modifications in the described device, and any further applications of the principles of the 
invention as described herein are contemplated as would normally occur to one skilled in the 
art to which the invention relates. 

10 FIG. 1 depicts system 1 10 of the present invention for fluorescent imaging of 

tissue 100. Tissue 100 has surface 101 and a heterogeneous composition as represented 
by regions 102, 103 underlying surface 101. Heterogeneities 102, 103 are generally 
not detectable by visual inspection of surface 101. 

System 1 1 0 includes modulated light source 120 to supply an intensity 

15 modulated excitation light of predetermined frequency and wavelength to tissue 100 
via optic fiber 123. Preferably, source 120 is a laser diode of conventional design with 
a modulated output in the 1 -500 MHz frequency range and a monochromatic output in 
the 100 to 1000 nm wavelength range. The specific wavelength is selected to excite a 
designated fluorophore in tissue 100. Beam splitter 126 may be employed to direct a 

20 small portion of the excitation signal to reference sensor 128 for processing purposes. 

System 1 10 also includes detection subsystem 140 which has optic fibers 143 to 
detect photons emitted from tissue 100 from a number of corresponding detection sites. 
Subsystem 140 includes one or more emission sensors 148. Detection subsystem 140 
also includes an interference filter to obtain a selected emission wavelength 

25 corresponding to emission of a designated fluorophore in tissue 100. In one 

embodiment, subsystem 140 includes a single sensor 148 and the signals from fibers 
143 are multiplexed. Preferably, sensors 128, 148 are Photo-multiplier Tubes (PMTs) 
or photodiodes but other sensor varieties, such as image intensifiers and charge-coupled 
devices, are also contemplated. 
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Sensors 128, 148 and source 120 are operatively coupled to heterodyne 
subsystem 130. Subsystem 130 is configured to obtain information about the phase, 
AC, and DC intensity of light detected with sensor 128 relative to light detected with 
the sensor 148 using conventional laser heterodyning techniques. In one embodiment, 
5 heterodyne subsystem 130 includes a signal synthesizer phase-locked to the repetition 
rate of a laser used for source 120. For this embodiment, subsystem 130 includes an 
amplifier to gain modulate sensors 128, 148 at a harmonic of a laser repetition rate 
(when a pulsed laser is used) or at the modulation frequency (when a modulated laser 
diode is used) plus an offset to provide the desired heterodyning. In one variation of 

10 this embodiment, an 80 MHz pulsed laser repetition rate is divided down to 1 0 MHz 
and input to the synthesizer, and a heterodyning offset of 1 00 kHz is input to the 
amplifiers for sensors 128, 148. 

Sensors 128, 148 are operatively coupled to processor 160. Processor 160 
includes input/control device 162, output device 164, and memory 166. Processor 160 

15 may be an electronic circuit comprised of one or more components. Similarly, 

processor 160 may be comprised of digital circuitry, analog circuitry, or both. Also, 
processor 1 60 may be programmable, an integrated state machine, or a hybrid 
combination thereof Preferably, input device 162 is a keyboard or input control of a 
conventional variety, and output device 166 is a Cathode Ray Tube (CRT) based video 

20 display, printer, or other image display system known to those skilled in the art. 
Memory 1 66 is preferably of the electronic (e.g. solid state), magnetic, or optical 
variety of the type readily available for use with electronic controllers or processors. 
Furthermore, Memory 166 may include an optical disk memory (CD), electromagnetic 
hard or floppy disk media, or a combination of these. 

25 FIG. 2 depicts one mode of operation of system 1 10 as process 210. Process 

210 includes mapping the spatial variation of fluorescence yield and lifetime with 
processor 1 60 and generating an image signal in accordance with the map. Output 
device 164 is configured to display an image in response to the image signal. Process 
210 begins by introducing a fluorescent contrast agent into tissue 100 in stage 212. 
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This agent provides a source of fluorescent emission for detection by subsystem 240. 
The configuration of the modulated light source 120, heterodyne subsystem 130, and 
detection subsystem 140 is designed to accommodate the excitation and emission 
properties of the selected fluorescent agent In other embodiments, endogenous 
fluorophores may be alternatively or additionally employed and system 1 10 adjusted 
accordingly. 

In stage 214, tissue 100 is excited by light source 120 configured according to 
the selected fluorophore. In stage 216, the phase, G obs , and log of AC intensity, M obs , 
of the emission at each detection site "i" relative to the excitation light from source 120 
are determined at the heterodyne (or offset) frequency. For "Di n number of detection 
sites, the detected or observed phase and AC intensity are indexed by "i" using the 
following notation: (e obs )j and (Mq^, respectively. Processor 1 60 stores the 
relative phase and AC intensity information in memory 166. 

In stage 21 8, a two dimensional grid is established for an area of tissue 100 
selected for imaging, and a matrix of grid points is established and indexed by "j". A 
uniform seed value for the fluorescent yield, yj = (Tm ax ^ m )j, and the fluorescent 
lifetime, (t)j, at each grid point j is assigned. These values are an initial homogeneous 
guess of the yield and lifetime values which are modified in later stages. The term 'V 
is the quantum efficiency of the fluorophore which varies with the environment of the 
surrounding of the fluorophore. The term " H ax ->m" is ^ ^sorption coefficient for 
the fluorophore and is the product of the extinction coefficient of the fluorophore based 
on the natural log and the concentration of the fluorophore. As a result, the yield, y = 
n^ax-nn' iS influenced by the surrounding metabolism and the uptake of the 
fluorophore. The uptake of certain known fluorophores vary with the type and 
condition of host tissue, providing another fluorescence characteristic useful to detect 
disease. The contrast provided by these properties is largely independent of 
fluorophore concentration. The initial estimate of fluorescent yield and lifetime are 
stored in memory 166 by processor 160 for later use. 
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After establishing this initial estimate of the fluorescence characteristics of 
yield, ^Hax^nr m ^ ^ et " ne ' T » processing loop 220 is entered in stage 230. 
Preferably, the stages of processing loop 220 are executed by processor 160 via 
preprogrammed software, dedicated hardware, or a combination of both as appropriate. 
5 To aid in understanding various mathematical aspects of -process 2 1 0 and loop 220, the 
following table of selected variables is listed: 



c velocity of light; 

D(r) optical diffusion coefficient; 

10 Di number of detection sites; 

f modulation frequency; 

I identity matrix; 

i detection site index; 

J Jacobian matrix relating the sensitivity at each grid point, j, to 

15 the response at each detection site; 

j grid point index; 

Jj i individual elements of the Jacobian matrix J; 

k source index; 

M log of AC intensity of modulated fluorescent light position; 

20 m index to multiple modulation frequencies; 

n average index of refraction; 

r position (in two or three dimensions); 

Sk number of modulated light sources; 

S(r, co) source term for the modulated light at position r and frequency 

25 co; 

Greek 

X 2 merit function representing the least squares error; 
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10 



15 



20 



<D x (r, co) 
Ha 

^a x -*c 
Ma x -> m 
M's 
9 

T 
CO 

Subscripts 

obs 

x 

m 



complex number representing photon flux in the frequency 
domain at position r and frequency co; 
quantum efficiency of fluorescent probe or dye; 
average absorption coefficient; 

absorption coefficient of the fluorescence light by both the non- 

fluorescing chromophores and fluorophore; 

absorption coefficient of the excitation light by both the non- 

fluorescing chromophores and fluorophore; 

adsorption coefficient due to non-fluorescing chromophores; 

adsorption coefficient of excitation light by fluorophores; 

effective scatting coefficient; 

phase-shift of one modulated light wave to another; 

lifetime of activated probe or dye at location r; 

angular modulation frequency, given by 2nf; 



observed or experimental data; 
excitation light; and 
fluorescence or emission light. 



In stage 230, phase and relative AC intensity at each detection site "i" is 
calculated as a function of the initial estimates of yield and lifetime for each grid point 
j. The calculated phase and intensity are represented at each detection site i as (8 m )j 
and (M m )j, respectively. The values for (6 m )j and (M m )j are determined using the 
25 diffusion equation approximation of the radiative transport equation. The diffusion 
equation approximation describes the spatial and temporal transport of light in tissues 
or multiply scattering media. A coupled frequency domain diffusion equation can be 
used to predict the excitation and emission fluence rates, <D x (r, o) and d> m (r, co), 
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respectively, at any location r within the selected grid of tissue 100 via equations (1) 
and (2): 

V • [D x (r)V<D x (r, ©)HHa x ( r )+i ©/Cnl^fr ©) +S x (r,co) = 0 

5 

V . P m (r)V(D m (r, <o)HHa m (r)+i <o/c n ]O m (r, co) +S m (r,o>) = 0 

The source term for the excitation light S x (r,co) is due to the sinusoidally modulated 
10 light at an angular frequency co = 2cof where f is typically in the MHz frequency range. 
The first term in both of the diffusion equations (1) and (2) represents the diffusive or 
"random-walk" transport of light where D x m is the optical diffusion coefficient of 
equation (3) as follows: 

15 D x>m =[3(M axm+ n' Sxm )]-l (3) . 

and |i a and |i' s are the absorption and isotropic scattering coefficients, respectively, for 
tissue 100, the medium of interest. The optical properties are dependent on the 
wavelength of light and thus are different for the excitation light from source 120 

20 (subscript x) and fluorescent emission detected with subsystem 140 (subscript m). The 
total absorption coefficient at the excitation wavelength, |i ax , is due to contributions 
from non-fluorescing chromophores as well as from fluorophores responsive to the 
excitation wavelength. The total absorption coefficient is given by the sum of 
absorption coefficients due to non-fluorescing chromophores, \i ax __± c , and fluorophores 

25 Ma x _»nr Generally it may be assumed that the absorption experienced at the 

fluorescent wavelength is due primarily to non-fluorescing chromophores. The velocity 
of light in tissue is Cn^dn where n is the average index of refraction. The source term 
for the fluorescent emission is dependent on the excitation light fluence, O x (r, co) and 
is given by equation (4) as follows: 
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S m M) - T 1 n ax ^ m (r)O x (r, co)[(l - icoT(r))/(l + ©2r(r)2)] (4) 

This term arises from the Fourier transform of the fluorescence decay term in the time 
5 domain following an incident pulse of excitation light where: t is the fluorophore 
lifetime, r\ is the quantum efficiency, and the absorption coefficient, Ha x -^m' * s * e 
product of the extinction coefficient based on natural log and the concentration of the 
fluorophore in the ground state. As previously indicated, the combined product, 
Wax^m' is termed * e fluorescent yield, y, and is proportional to the generated 

10 fluorescence fluence. Substitution of equation (4) into equation (2) facilitates 

determination of <D m for each grid point "j." The solution of the diffusion equations 
(1) and (2) for the two-dimensional area defined by the grid points "j" may be readily 
extended to three dimensions to estimate spatial variation of one or more fluorescence 
characteristics in a selected volume with "r" corresponding to position in three 

15 dimensions. 

Both diffusion equations (1) and (2) are linear complex elliptic equations that 
can be solved as boundary value problems for the complex quantities <D x (r, co) and 
<l> m (r, co). This solution employs the method of finite differences to create 
corresponding finite difference equations. These difference equations are utilized to 

20 obtain an approximate solution at each grid point, j. This method of solution is 

described in other contexts in Fulton et al., Multigrid Method for Elli ptic Problems. A 
Review . 1 14 American Meteorological Society pp. 943-59 (May 1986); and B.W. 
PogueetaL, Initial Assessment of a Simple System for Freque ncy Domain Diffuse 
Optical Tomography . 40 Physics in Medicine and Biology pp. 1709- 1729 (1995). One 

25 preferred method of performing this solution is with the MUDPACK routines described 
in Adams, J.C., MIJDPACK; Multigrid Portable Fortran Software for the Efficient 
Solution of Linear Ellip tic Partial Differential Equations . 34 App. Math Comp. p. 133 
(1989). For the solution of the diffusion equations, it is assumed that ^m^xfo 05 ) = 0 on 
the surface 101 of tissue 100 which is known as the zero fluence boundary condition. It 
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should be recognized that other boundary conditions may be selected and the method of 
solution varied accordingly. 

The diffusion equations (1) and (2) may be solved for a complex number for 
G> m at each grid point, j. The detected signal at the surface is proportional to the 
5 normal component of the gradient of the photon fluence. To approximate the signal at 
detector site "i" located on surface 101 of tissue 100, the d> m value at an internal grid 
point closest to the site is selected which follows from the relationship that the normal 
component of the photon fluence gradient is proportional to <t> m just inside surface 101. 
The calculated phase-lag, 0 m , and the log of AC intensity, M m , at the detection sites 

10 "Di" are calculated from the imaginary and real parts of the complex <E> m value with 
respect to the phase and the AC intensity of source 120. 

The diffusion equations (1) and (2) provide insight into the sensitivity of 
changing the fluorescent optical properties of tissue 100 on 0 m and M m measured at 
the detector sites i. This insight results from a series of calculations while fixing 

15 various parameters of the diffusion equations (1) and (2). These calculations assume 
circular tissue phantom 300 with an embedded, heterogeneity 302 hidden in phantom 
background 303 as illustrated in FIG. 3. A two-dimensional grid is established for 
phantom 300 and may easily be expanded to three dimensions. Under these simulated 
conditions, a large value is assigned to absorption coefficients for both excitation and 

20 fluorescent light at all grid points outside the simulated tissue phantom. The four 
sources S1-S4 of FIG. 3 (Sk = 4) are simulated by assigning an arbitrary complex 
number at a grid point near the surface closest to each source. The twenty detection 
sites Dl -D20 of FIG. 3 (Di = 20) are simulated by using the calculated values 
determined from <D m at the grid point "j" closest to the detection site. The simulated 

25 solutions to diffusion equations (1) and (2) were obtained in two dimensions for a 65 x 
65 grid covering a 100mm diameter circular tissue phantom 300 with a circular, 
embedded heterogeneity having a 30mm diameter and located at the center of the tissue 
phantom 300 (this location differs slightly from the configuration of heterogeneity 302 
of FIG, 3). The simulated measurements of fluorescent phase-shift and AC intensity 
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10 



are reported for 20, equally spaced, circumferentially located detection sites D1-D20. 
The modulation frequency, f, was set equal to 150 MHz. The optical properties of the 
heterogeneity and the background are shown in Table 1 as follows: 







^x or 






X 


frequency 


(mm" 1 ) 


(mm* 1 ) 




(mm* 1 ) 


background 


background 


(MHz) 






(mm* 1 ) 




(mm* 1 ) 


(ns) 






0.0 


1.0 


0.0 


1.0 x 10" 5 


1.0 


150.0 

















Table 1 



In order to evaluate the influence of Wax— ►m' ® m ^ m were com P ute ^ at 
each detection site D1-D20 as the value of "lHax-Mn "* ^ e heterogeneity increased 
from 10" 4 mm"' to 10"' mm"' and as "IM-ax— >m * n t ^ ie background 303 was maintained 

15 constant. The lifetime, x, was set equal to 1 ns for both the object and the background 
causing contrast due to differences in "nii ax ^ m . The plots of 6 m and M m are shown in 
FIGS. 4 and 5 respectively for one active source SI. As ^n^ax-^m °f heterogeneity 102 
increases to higher values, the AC intensity approaches an upper limit similar to what is 
expected in dilute non-scattering solutions. FIG. 5 shows how the fluorescent phase- 

20 shift, 0 m , decreases as the absorption coefficient due to the fluorophore, Ha x ->m IS 
decreased 10 to 100 times the background. From these simulations, M m appears to be 
directly dependent upon changes in ilHax-nn °^ a s * mu ' ate d tissue heterogeneity 102 
whereas G m is indirectly dependent on TiHa x -*m ^ ue t0 c ' ian S es * n photon migration. 
In order to evaluate the influence of x, 8 m and M m were calculated at each 

25 detection site D1-D20 as the values of x in the heterogeneity varied from 10"' ns to 10^ 
ns and the value of x in the background was held at 1 ns. The background *lHax-Mn 
was set to 10'^mm"' and *lHax-»m for the heterogeneity was set to 10"^mm"'. As 
shown in Fig. 6, the detected AC intensity increases as x decreases. Fig. 7 illustrates 
the values of the fluorescent phase-shift at.each detection site as the lifetime of the 
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heterogeneity is changed from 0.1 ns to 1000 ns. At a given modulation frequency (150 
MHz in this calculation), 0 m first increases, reaches a maximum and then subsequently 
decreases as t is increased from 0. 1 ns to 1000 ns. Therefore, both 9 m and M m at each 
detection site D1-D20 appear to be directly influenced by the value of lifetime in the 
heterogeneity. 

Referring back to FIG. 2, in stage 240, the calculated emission phase and 
intensity, (Q m )i and (M m )j, are compared to the measured emission phase and intensity, 
( e obs)i ^ (M 0 b s )i, for each detection site "i" to identify a difference or "error" 
between the measured and calculated values. Because OlHax-Hn^j im P acts ( M m)i » ^ s 
comparison is posed in the form of the merit function %y} of equation (5) as follows: 

Sk Di 

X M 2 = (1/Sk) I (1/Di) E [((Mob^.CM^iyoM^ ( 5 ) 
k=l i=l 

where aj^ is the typical standard deviation of noise in M m , taken to be 0.01 ; Sk = 
number of excitation source sites indexed by k; and Di = number of detection sites 
indexed by i. The goal of the algorithm is to minimize %^ by appropriate updates of 
(^ax-* m^j* A ^ er 311 i 1 "^ update of (Tm ax-> m )j another merit function in terms of 
(t)j participates in the comparison of stage 240. This merit function, Xx^ » * s presented 
as equation (6) as follows: 

Sk Di 

X T 2 = (1/Sk) I (1/Di) Z [((Mobs)i-(^ (6) 
k=l i=l 

where cq is the typical standard deviation of noise in (8 m )j, taken to be 1 degree; Sk = 
number of excitation source sites indexed to k; and Di = number of detection sites 
indexed to i. Since the lifetime influences both (0 m )i and (M m )j, the phase and AC 
intensity are used in equation (6). 
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After the comparison of stage 240 is performed by calculating the merit 
functions > control flows to conditional 250 to test whether the comparison of 

the observed values, (9 0 bs)i (Mobs)i> t0 * e calculated values (9 m )i and (M m )j 
via the merit functions meets a selected convergence criteria. This criteria corresponds 
5 to the degree of tolerable error in determining the yield and lifetime values. For one 
embodiment, convergence is achieved when any of the following three quantities, (i) 
X^, (ii) change in $ in successive iterations of loop 220, (ii) relative change in yp- in 
successive iterations of loop 220 is less than a predetermined threshold value of 1 .0 x 
1 0"^. In other embodiments a different comparison calculation and associated 

10 conditional may be employed as would occur to one skilled in the art. If conditional 
250 is satisfied, control flows to stage 270 and loop 220 is exited; however, if the 
criteria is not satisfied, execution of loop 220 continues in stage 260. 

In stage 260, the yield, (y)j = CnHax-^m^p ^ lifetime, (x)j, for each grid point 
j is updated so that these values may reach the minimum error corresponding to the 

15 comparison stage 240 and conditional 250 test. In order to update these values, * 
Jacobian matrices are used which describe the sensitivity of the response at each 
detection position i to changes in (y)j = Omax-nn^j' ^ ^fe^ 1116 ' ( T )j> at eac ^ E"d 
point, j. Three Jacobian matrices are employed: J (M, *lHax-» 311(1 
t). The elements Ji j of these Jacobian matricies are given by Ji j = 

20 [5Mj/(a(Tm ax ^ m )j]; JiJ = [dMj/dtj]; and Ji j = [dQ/dxj], respectively. These elements 
may be calculated by solving the diffusion (1) and (2) four times for each grid point, j 
to obtain M m ^ and G m ,j calculated with (x )j and (t+ 5t)j and with OlHax-nn^j ^ 
Omax-nn + ^^ax-^m^T From ' east squares minimization, the update to yield and 
lifetime is calculated. In one preferred embodiment, this updating algorithm is adapted 
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from an algorithm used to reconstruct images obtained by electrical impedance 
tomography like the algorithm suggested by Yorkey, et aL, Comparing reconstruction 
Algorithms for Electrical Impedance Tomographv T 34 Transactions in Biomedical 
Engineering pp. 843-52 (1987). The Jacobian Matrices are used to solve for update 
vectors, [Arm ax _ >m ] and [Ax], to estimated yield and lifetime vectors, [r\\Xa x _+ m ] 
and [t], respectively. These vectors are of a dimension corresponding to the number of 
grid points. At each iteration through loop 220, the following Jacobian equations (7) 
and (8) are solved to determine the update for the estimated yield and lifetime vectors: 



J(M.Tm ) T JfM.nu ) 



x "» m +31 1 



J M 



Ka_J = 



x — > m 



■(Mm^ -M m ) 



M 



(7) 



J(M.T) T J(M.T) J(6,T) T J(e,T) 



M 



R- 



j(m.t) t — - J(e.t) T - 

■IMm . -Mm) + r-(e m , -e m ) 



M 



obs 



obs 



(8) 



M m ^ bs and M m are the observed and calculated vectors of the log of AC intensity at 
each of the i detection sites, respectively, and 9 m are the observed and 
calculated vectors of the phase lag at each of the i detection sites, respectively . Due to 
the ill-conditioned nature of the Jacobian matrices, the terms X\l or 'hqi are added as 
part of a Marquardt minimization scheme where I is an identity matrix. The parameters 
or A.2 are adjusted via a Maquardt-Levenberg type algorithm of the type disclosed in 
Press et al., Numerical Recipes: The Art of Scien tific Computing . (Cambridge 
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University Press, 1992). Conventional numerical methods are employed to solve the 
simultaneous linear algebraic equations resulting from the Jacobian matrix equations 
(7) and (8). The Jacobian matrices are re-calculated at each iteration through loop 220. 
It has been found that equations (7) and (8) provide a way to select appropriate changes 
5 to the yield and lifetime estimates; however, other numerical approaches to recursively 
iterate to acceptable estimates as would occur to one skilled in the art are also 
contemplated. Once the update is complete, control returns to stage 230. 

If the convergence criteria is satisfied in conditional 250, then estimation of 
yield and lifetime for the grid points has reached an acceptable minimum and control 
10 flows to stage 270. In stage 270 an image signal is generated by processor 1 60 from the 
spatial variation of the yield and/or lifetime fluorescence characteristics. This image 
signal is sent to output device 164 which displays an image in response. Because the 
fluorescence characteristics of yield and lifetime typically vary with the biologic 
environment of the fluorophore, this image is generally indicative of tissue variation 
15 and offers the capability to detect heterogeneities 102, 103. For example, laser diodes 
capable of supplying Near infrared (NIR) light that can penetrate tissue several 
centimeters, and fluorescent contrast agents responsive to NIR light may be used to 
provide a viable imaging system. In one embodiment, this system is adapted for use 
with an endoscope. 

20 Besides yield and lifetime, the spatial variation of other fluorescence 

characteristics useful to distinguish diseased tissues may be mapped using the diffusion 
equations (1) and (2). Such alternative fluorescence characteristics include, but are not 
limited to, quantum efficiency r\ and/or fluorescent absorption coefficient Hax-»m 
determined as separate properties independent of the yield product. 

25 In another embodiment of the present invention, the photon fluence equation 

and Jacobian estimation process is adapted to determine a map of a designated 
fluorophore uptake concentration. For this embodiment, a first map of chromophore 
adsorption coefficients Ma x -)-c scarring coefficients are determined in the 
absence of the designated fluorophore by estimating the chromophore adsorption 
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coefficient H ax ->c ^ scattering coefl5cient ^i' s at each grid point j in place of the 
yield and lifetime estimates. Diffusion equation (1) for <D x (r, may be used in 
conjunction with modified Jacobian equations (7) and (8) to create this first map. The 
modification substitutes the chromophore adsorption and scattering coefficients in place 
5 of the yield and after adaptation to accommodate these new characteristics as follows: 



10 




and 



15 



20 



J<M x ,n s ) T - 



(Mx obs -M X ), + 



J 8 



-(ex obs -e x ) 



(10) 



25 



The elements of the four Jacobian matrices employed, J(M x ,n a ), J(M x ,m), J^XiH ). 

respectively. Updates to the absorption and scattering map were conducted to minimize the 
merit function x^: 



n s k = 1 n d i=l 



xobs,i m x,j 



xobs,i~ e x,j 



(11) 



where n, — Sk and = Di. 
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After generating the first map, the designated fluorescent contrast agent is 
introduced, and the total adsorption coefficient ji ax is determined by substituting ji ax 
in place of H ax -*c m ^^ons (9)-(l 1) to obtain a second map of the total adsorption 
coefficient. Noting that ji ax = H ax _» m + n ax-> c , and that the uptake of the 
fluorescing contrast agent is directly proportional to H ax _> m , u P take concentration may 
be mapped by determining a difference between the adsorption coefficient variations 
for the first and second maps. This "difference map" may then be used to generate an 
image corresponding to the uptake concentration. 

Another alternative embodiment measures the emission responsive to each of a 
number of light source modulation frequencies f. The total number of different 
frequencies employed is designated Mf. To obtain this additional data, an iteration of 
loop 220 is performed for each frequency f indexed to m. The number of sources, Sk 
and detection sites Di are indexed to k and i, respectively. This additional data may be 
used to enhance imaging results obtained with system 1 10 or to permit reduction of the 
number of detection sites or excitation source sites in the evaluation. A representative 
merit function corresponding to this additional data is given in equation (12) as follows: 

Xr 2 = 

Mf Sk Di 

(1/Mf) Z (1/Sk) Z (1/Di) Z [((M obs )i - (M m )i)/a M ]2 + [((Q ohs)i - (e m )i)/a e ]2 (12) 
m=l k=l i=l 

Besides fluorescence yield and lifetime, the multi-frequency method can be employed 
to map other optical characteristics of interest. Besides a sinusoidally modulated light 
source, the present invention may be adapted to operate with a pulsed or other time- 
varying excitation light source in alternative embodiments. 

FIG. 15 depicts an optical system 410 of another embodiment of the present 
invention. This system includes modulated light source 420 with laser driver 422, 
operatively coupled laser diode 424, and reference frequency generator 426. Source 
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420 is configured to deliver modulated light to tissue phantom 400, and the re-emitted 
light from the phantom is focused onto a gain modulated image intensifier 430 via 
50mm lens 432. Intensifier 430 includes a photocathode face which converts photons 
to electrons, a Multi-Channel Plate (MCP) which multiplies the electronic signal by 
5 avalanche multiplication, and a phosphorescent screen which converts electrons into an 
optical image. Preferably, intensifier 430 is a fast intensifier, of the variety 
manufactured by Litton Electronics, Inc., which enables modulation by applying a DC 
bias and an RF signal from amplifier 428 between the photocathode and the MCP. For 
this example, the modulation of the image from intensifier 430 is phase-locked to the 
10 laser diode 424 by a 10 MHz output signal from synthesizer 426. By modulating the 
laser diode 424 and the image intensifier 430 at the same frequency, a steady-state 
image results on the phosphor screen. U.S. Patent No. 5,213,105 to Gratton et al. 
provides additional background concerning certain aspects of this technique. The 
image from the phosphor screen is focused through interference filter 433 on a Charge 
15 Coupled Device (CCD) camera 434 via 150 mm macro lens 436. Camera 434 has a 
512x512 array of CCD detectors configured to provide a corresponding pixelated 
image. Camera 434 is operatively coupled to processor 460 of a similar configuration 
to processor 160 previously described. 

Following each acquired image, a phase delay between the image intensifier 430 
20 and the laser diode 424 is induced by stepping the phase of the image intensifier 430 to 
values between 0 and 360 degrees with the frequency synthesizer 452 under the control 
of processor 460. Since the gain modulation of image intensifier 430 and laser diode 
424 occurs at the same frequency, homodyning results in a steady phosphorescent 
image on intensifier 430 which is dependent upon phase. Preferably, control between 
25 synthesizer 452 and processor 460 is obtained by a conventional GPIB interface. 

Images from the phosphorescent screen of the image intensifier 430 are then gathered at 
each phase delay. The incremental phase delayed images are then used to generate a 
map of phase-shift and intensity modulation ratio between the excitation and re-emitted 
light from phantom 400. By applying interference or appropriate optical 
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filters, the emission light may be selectively separated from the excitation light and 
measured. Camera 434 output may be processed by processor 460 using process 210. 

The present invention will be further described with reference to the following 
specific examples 1-3. It will be understood that these examples are illustrative and not 
restrictive in nature. Examples 1-3 involve the computer simulation of the process 
210. Simulations of this kind, including the simulation of tissue, are an acceptable 
means of demonstrating fluorescent spectroscopic imaging performance to those skilled 
in the art. The examples use simulated values obtained by solving the diffusion 
equations (1) and (2 ) for G m and M m under the conditions of table 2 as follows: 



Case 






^ or 


X 


ra i-im 


Gaussian 


Gaussian 










(background) 


(background) 


Noise 
in log of AC 


Noise 
in phase 




(mm* 1 ) 


(nun -1 ) 


(mm-') 


(ns) 


(mm -1 ) 


intensity 
<?M 


ce (degrees) 


5.1 


0.0 


0.0 


1.0 


10.0 


1.0 x 10-5 


0.01 


0.1 


5.2 


1.0 x 10-3 


0.0 


1.0 


10.0 


1.0 x 10-5 


0.01 


0.1 


5.3 


0.0 


0.0 


1.0 


10.0 


1.0 x 10- 5 


0.01 


1.0 



Table 2 

The examples simulate tissue phantom 300 of FIG. 3 having a 100mm diameter. 
Values of 9 m and M m were computed at each of the D1-D20 detection sites of FIG. 3 
in response to the 4 modulated light sources S1-S4 located at the periphery. The 
excitation light modulation frequency f was simulated at 150 MHz. Diffusion 
equations (1) and (2) were solved to provide 80 simulated values of 6 m and M m 
corresponding to the various combinations of detection and source sites (Sk * Di = 
4x20 = 80). Gaussian noise with a standard deviation of 0.1 degrees (or a liberal 1 
degree) in 9 m and 1% in M m were superimposed on the diffusion equation solutions. 
Adapted MUDPACK routines were used to solve the diffusion equations (1) and (2) on 
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a SunSparclO computer. These obtained data sets were used as simulated input data to 
process 210 for examples 1-3. The results are shown in tables 3 and 4 are as follows: 



5 


Case 


Area, object 1 
(mm 2 ) 


Location, object 1 
(x f y), (mm, mm) 


Area, object 2 
(mm 2 ) 


Location, object 2 
(x,y), (mm, mm) 




5.1 


706.0 (expected) 
742.2 (obtained) 


(60,60) (expected) 
(60.8,58.5) (obtained) 


not applicable 


not applicable 


10 


5.2 


706.0 (expected) 

703.1 (obtained) 


(60,60) (expected) 
(59.4,58.3) (obtained) 


not applicable 


not applicable 




5.3 


314.1 (expected) 
381.0 (obtained) 


(32.3,67.7) (expected) 
(34.0,67.7) (obtained) 


314.1 (expected) 
342.0 (obtained) 


(67.7,32.3) (expected) 
(65.0,35.0) (obtained) 



Table 3 

15 



25 



Case 


(object) 


x (object) 




(mm* 1 ) 


(ns) 


5.1 


1.0 x 10' 3 (expected) 


1 .0 (expected) 




0.93 x 10* 3 (obtained) 


1.03 (obtained) 


5.2 


1.0 x 10* 3 (expected) 


1.0 (expected) 




0.8 x 10* 3 (obtained) 


0,7 (obtained) 


5.3 


(top left object): 


(top left object): 




1.0 x 10- 3 (expected) 


1.0 (expected) 




2 x 10* 3 (obtained) 


4.1 (obtained) 




(bottom right object): 


(bottom right object): 




2.0 x 10- 3 (expected) 


2.0 (expected) 




1.8 x 10- 3 (obtained) 


3.5 (obtained) 



Table 4 
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Example 1 

Example 1 reconstructs fluorescent yield and lifetime with no absorption due to 
non-fluorescing chromophores. To simulate the experimental data for this example, 
the fluorescent yield, OlHax-Wj' for ^ background and the heterogeneity 302 were 
chosen as 1 x 10' 5 mm" 1 and 1 x 10" 3 mm" 1 respectively and the fluorescence lifetime, 
(x )j, for the background and the heterogeneity 302 chosen as 10 ns and 1 ns 
respectively. During the execution of loop 220, no a priori knowledge of either the 
heterogeneity 302 location or the background fluorescence properties was assumed and 
a uniform guess of Ix 10* 5 mm' 1 and 10 ns was given for the fluorescence yield, 
(Tl^ax-^m^j' md Iifet »me, ( T )y respectively. Convergence was achieved in less than 50 
iterations of Loop 220 (computational time on a SunSparclO: 2 hours) for a two 
dimensional 17x17 grid. The average values of Tm ax -Mn and x in the grid points 
which occupy the simulated object converge within 50 iterations to r\\i ax _+ = 0.93 x 
10° mm" 1 and t =1.03 ns are illustrated in FIGS. 8 and 9, respectively. FIGS. 1 0 and 
1 1 illustrate the reconstructed images from the mapped values of Tin ax _ >m and t, 
respectively, and are representative of the expected images. The images were smoothed 
by interpolation in examples 1-3 to remove spurious points which had unphysically 
high values, but were surround by values within a physically achievable range. These 
spurious values were replaced by the average background fluorescence yield and 
lifetime obtained from simulation of loop 220. 

The average values of im ax _^ m in the grid points which occupy the simulated 
background converge within 50 iterations to 9 x 10' 5 mm' 1 . The value of the 
background converges to 5.4 ns. The dependence of the final images on the choice of 
the initial guess was examined by providing an initial uniform guess of 1 x 10" 4 mm' 1 
and 10ns for (^ ax ^ m )j, md lifetime > ( T )j> respectively. This resulted in similar 
images to those obtained in FIGS, 1 0 and 1 1 . 

The location of heterogeneity 302 was identified as consisting of all the grid 
points with TlH ax ^ m higher than 35% (arbitrarily chosen) of the peak value of the 
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TlHax-mi (F" 10 - 10 )- The average of the coordinates of all the identified object grid 
points was the position (60.8, 58.5) which is close to position (60, 60) that was used to 
simulate the experimental data. As listed in Table 3, the area of the heterogeneity based 
upon our arbitrary definition for identification was 72 mm 2 , close to that used to 
generate our simulated experimental data. 



Example 2 

Example 2 reconstructs fluorescent yield and lifetime with a simulated 
chromophore absorption configured to mimic tissue. The same hidden heterogeneity as 
well as optical parameters and simulation equipment were used as described in 
Example 1 except that a uniform background chromophore absorption coefficient, 
^a x -> of 1 xl0 ~ 3 mm " 1 was used to generate the simulated experimental data. While 
excitation light propagation was not employed for image reconstruction, we considered 
this optical property known to estimate the best possible performance for inverse image 
reconstruction under physiological conditions. The two-dimensional reconstructed 
spatial map of the fluorescence yield, (Tm ax -> m )j' 211(1 Iifetin ie, (x)j, are shown in 
FIGS. 12 and 13, respectively. As shown in Table 3, the mean value of location of the 
object according to our criterion based on Tm ax ->m occurred position (59.4, 58.3) 
consistent with the conditions used to simulate the experimental data. The dimension 
of the heterogeneity based upon our arbitrary definition for identification (all grid 
points with TiHa x -> m tegfer than 35% of the maximum) were 703 mm 2 which is close 
to that used to generate our simulated experimental data. The average values of 
Wax-nn 311(1 T in ^ S rid P oints which occupy the simulated object converge within 
50 iterations to the values of r|jx ax _ >m = 0.8 xlO" 3 mm -1 and x = 0.7 ns consistent 
with the values used to generate the simulated experimental data (see Table 3). The 
average values of i1Hax->m ^ T in ^ S rid P° ints which occupy the simulated 
background converge within 50 iterations to values similar to that reported for Example 
1. 
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Example 3 

Example 3 simulated two hidden heterogeneities in the tissue phantom (not 
5 shown in FIG. 3). In this case, the same optical parameters were used as described in 
example 1 except that the fluorescence yield Tm ax _^ m for the objects 1 and 2 was 
chosen as 1 x 10"3mm'l and 2 x 10'^mm'^ respectively and lifetime x for the 
heterogeneities chosen as 1 ns and 2 ns, respectively. A 33 x33 grid was employed 
instead of a 17 xl7 grid. An image corresponding to the mapping of yield is depicted 
10 in FIG. 14. 

All publications and patent applications cited in this specification are herein 
incorporated by reference as if each individual publication or patent application were 
specifically and individually indicated to be incorporated by reference. While the 
15 invention has been illustrated and described in detail in the drawings and foregoing 

description, the same is to be considered as illustrative and not restrictive in character, it 
being understood that only the preferred embodiment has been shown and described 
and that all changes and modifications that come within the spirit of the invention are 
desired to be protected. 
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CLAIMS 

What is claimed is: 

1 . A method of imaging, comprising the steps of: 

5 (a) exposing a surface of a light scattering material to an excitation light from a light 

source, the material having a heterogeneous composition under the surface; 

(b) detecting a fluorescent emission from the material in response to step (a); 

(c) establishing an estimate of spatial variation of a fluorescence characteristic of the 
material; 

10 (d) determining a calculated emission as a function of the estimate; 

(e) comparing the calculated emission to the emission detected in step (b) to determine 
an error; 

(£) providing a modified estimate of spatial variation of the fluorescence characteristic 
and repeating steps (d) through (f) until the error reaches a desired minimum; and 
15 (g) generating an image of the material from the modified estimate, the image 

corresponding to the heterogeneous composition of the material. 

2. The method of claim 1, further comprising introducing a fluorescent contrast agent into 
the material. 

20 

3. The method of claim 1 , wherein step (f) includes applying a Jacobian matrix. 

4. The method of claim 1 , wherein the expected emission is determined as a function of a 
diffusion equation, 

25 

5. The method of claim 1, wherein step (ej includes comparing intensity and phase of the 
calculated emission to intensity and phase of the emission detected in step (b). 
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6. The method of claim 1, wherein the material includes a fluorescent contrast agent and 
the fluorescence characteristic is a function of at least one of fluorescence quantum 
efficiency, fluorescence lifetime, and concentration of the fluorescent contrast agent. 

5 7. The method of claim 6, wherein the fluorescence characteristic is a function of 

fluorescence quantum efficiency, the light source is intensity modulated at a predetermined 
frequency, step (e) includes comparing AC intensity and phase of the calculated emission to 
intensity and phase of the emission detected in step (b), step (f) includes applying a Jacobian 
operator, and the expected emission is determined as a function of a photon fluence rate 
10 within the tissue. 

8. A method of imaging, comprising: 

(a) exposing a living biologic tissue to an excitation light from a light source; 

(b) detecting a fluorescent emission from the tissue in response to said exposing; 

15 ( c ) determining a spatial variation of a fluorescence characteristic within the tissue as a 

function of the emission with a processor; and 

(d) generating an image of the tissue in accordance with the spatial variation. 

9. The method of claim 8, further comprising introducing a fluorophore into the tissue. 

20 

10. The method of claim 8, wherein the fluorescence characteristic is a function of at least 
one of fluorescence lifetime, fluorescence quantum efficiency, or fluorescence absorption. 

1 1 . The method of claim 9, wherein said exposing includes positioning the source adjacent 
25 a surface of the tissue, and said detecting includes sensing emitted light at the surface. 

12. The method of claim 8, wherein said exposing includes positioning a number of 
modulated light sources adjacent a surface of the tissue and said detecting includes sensing 
emitted light at a number of locations along the surface. 
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13. The method of claim 8, wherein said determining includes (i) establishing an estimate 
of spatial variation of the fluorescence characteristic, (ii) determining a calculated emission as 
a function of the estimate, (iii) comparing the calculated emission to the fluorescent emission 
of said detecting to determine an error, (iv) providing a modified estimate of spatial variation 
of the fluorescence characteristic, and (v) repeating said determining and said comparing with 
the modified estimate until the error reaches a desired minimum. 

1 4. A method of imaging, comprising: 

(a) introducing a fluorescent agent into a light scattering biologic tissue with a 
heterogeneous composition; 

(b) exposing a surface of the tissue to light from a light source to excite the agent; 

(c) detecting an emission of light from the tissue in response to said exposing; 

(d) determining a number of values representative of a fluorescence characteristic of 
the tissue as a function of position and the emission; and 

(e) generating an image in accordance with the values, the image corresponding to the 
heterogeneous composition of the tissue. 

15. The method of claim 14, wherein the fluorescence characteristic is a function of 
fluorescence quantum efficiency. 

1 6. The method of claim 1 4, wherein the fluorescence characteristic corresponds to 
fluorescence lifetime. 

17. The method of claim 14, wherein said detecting includes sensing fluorescent light 
emitted in response to said exposing at a number of locations along the surface. 

18. The method of claim 14, wherein said exposing includes providing light from a number 
of sources spaced along the surface. 
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19. The method of claim 14, wherein said exposing includes exciting the fluorescent agent 
in the tissue with modulated light at a number of different frequencies, and the values are 
determined as a function of the different frequencies. 

5 

20. The method of claim 14, wherein said determining includes (i) establishing an estimate 
of the values, (ii) determining a calculated emission as a function of the estimate, (iii) 
comparing the calculated emission to the emission of said detecting to determine an error, (iv) 
providing a modified estimate of the values, and (v) repeating said determining and said 

10 comparing with the modified estimate until the error reaches a desired minimum. 

21. A system for imaging a light scattering tissue having a heterogeneous composition and 
containing a fluorophore, comprising: 

(a) a light source adapted to excite the fluorophore; 
15 (b) a sensor configured to provide a detected light signal corresponding to a fluorescent 

emission from the tissue in response to light from said source; 

(c) a processor operatively coupled to said sensor and responsive to said detected light 
signal to provide a number of values representative of a fluorescence characteristic of the 
tissue as a function of position, the fluorescence characteristic corresponding to at least one of 

20 fluorescence lifetime, fluorescence quantum efficiency, and fluorescence absorption, said 
processor being configured to generate an image signal as a function of said values; and 

(d) an output device responsive to said image signal to provide an image corresponding 
to the heterogeneous composition of the tissue. 

25 22. The system of claim 2 1 , further comprising a number of modulated light sources. 

23. The system of claim 2 1 wherein said fluorescence characteristic corresponds to 
fluorescent yield. 
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24. The system of claim 21, wherein said sensor is configured to detect said emission at a 
number of locations along a surface of the tissue. 

25. The system of claim 21, wherein said processor is configured to determine said values 
from a comparison of a calculated emission to an observed emission derived from said 
detected light signal, said calculated emission being determined as a function of an estimated 
spatial variation of said fluorescence characteristic, said estimated variation being updated 
and said comparison being repeated until a difference between said calculated emission and 
said observed emission reaches a desired minimum. 

26. The system of claim 2 1 , wherein said source includes a laser diode and said sensor 
includes a CCD camera. 
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